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The next-to- leading term of the renormalized stress-energy tensor of the quantized massive 

field with an arbitrary curvature coupling in the spacetime of the Schwarzschild black hole 

is constructed. It is achieved by functional differentiation of the DeWitt-Schwinger effective 
O 

action involving coincidence limit of the Hadamard- Mmakshisundaram-DeWitt-Seely coef- 



ficients 03 and 04. The back reaction of the quantized field upon the Schwarzschild black 
hole is briefly discussed. 



PACS numbers: 04.62. +v, 04.70.-s 



INTRODUCTION 



If the Compton length, A c = K/mc, associated with a quantized massive field is much smaller 



than a characteristic radius of curvature, L, (where the latter means, as usual, any length scale 
of the background geometry) then the nonlocal contribution to the renormalized effective action, 
Wr, can be neglected and its series expansion in m~ 2 can be constructed using DeWitt-Schwinger 
method. Since in the renormalization prescription one has to absorb the first three terms of the 
expansion into the classical action of the quadratic gravity with the cosmological term, the lowest 
nonvanishing term of the Wr is to be constructed from the (integrated) coincidence limit of the 
fourth Hadamard- Minakshisundaram-DeWitt-Seely coefficient, [0,3], whereas next to leading term 
is constructed form [04]. Generally one has 



W " 32, 

n=3 



htfl (^2)n- 2 / ^ x VaK]. (i) 
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For the technical details of this approach the reader is referred, for example, to Refs. and the 
references cited therein. 

It is a well known fact that the complexity of [a n ] increases rapidly with n making calculations 
of the coefficients for n > 2 a highly nontrivial task. It is expected therefore that the applicability 
of the series ([I]), truncated at some definite n, will be limited to the simplest geometries with 
symmetries. On the other hand, however, as the coefficients depend on the background geometry, 
and, possibly, on a "potential" term, they can be used to construct the renormalized stress-energy 
tensor, T%, by functional differentiation of Wr with respect to the metric. Such a tensor can be 
defined in a wide class of geometries, and, by construction, it gives a unique opportunity to study 
the back reaction on the metric in a self-consistent way. Of course, the results of such calculations 
should be interpreted with care as the particle creation, which is a nonlocal process, is ignored. 

The coefficient [aol has been calculated by DeWitt 131 whereas [as] by has been obtained by Sakai 

n n n □ q 

and Gilkey [J, |5j; the fifth coefficient, [04], has been calculated in Refs. [g, |7], |8j|. The results for 
[04] are rather hard to compare as there are various simplification strategies that can be employed, 
and, unfortunately, some of the results contain not only typographical errors. Moreover, a compact 
or even tricky notation is of little help in situations when the main task is to calculate the stress- 
energy tensor in a specific spacetime. Therefore, in order to construct the approximation to the 
renormalized stress-energy tensor we have independently calculated [04] for a massive scalar field 
with an arbitrary curvature coupling satisfying the equation 

(-□ + £R + m 2 )0 = O, (2) 

where £ is the dimensionless parameter describing the curvature coupling and R is the curvature 



scalar, using fully covariant method of DeWitt 



the Riemann normal coordinates 



and checked the calculations constructing [04] in 



9]. The thus calculated coefficients have been compared among 



themselves and with their known values in concrete geometries. For example, when specialized 
to n = 4, the coefficient [04] precisely reproduces the coefficient obtained from Dowker's general 
formula for [a n ] in the de Sitter spacetime [l(J 

dS 6 ±\(2^-l)B 2k \ 



^ S ( 4 « 2 ) 4 ^ fc!(4-fc)! 



(3) 



105a 8 ' 

where B 2 h are Bernoulli numbers and a is the radius of the curvature. It is zero in the optical 
version of the Nariai metric, as expected. Moreover, as an additional partial check, we have also 
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calculated the basic ingredient of the DeWitt method, [n 5 c], in two different ways, where the 
biscalar a(x,x') is half the square of the geodetic distance between points x and x' . Subsequently, 
making use of the standard formula 

7- = ±™« (4) 
og a b 

we have constructed the next-to-leading (i.e. m -4 ) term of the renormalized stress-energy tensor 
in a general spacetime. To the best of our knowledge it is the first attempt to go beyond the first 
order (i.e. m~ 2 ) in the calculations of this type. 

The DeWitt method is easily programmable, and the number of terms that appear at intermedi- 
ate stages of calculations can be reduced significantly by a carefully chosen simplification strategy. 
On the other hand, the calculations carried out in the Riemann normal coordinates are extremely 



fast 33]. The calculations of the coefficient [04] and its functional derivatives with respect to 
the metric tensor have been carried out with the aid of FORM 111] and its multithread version 
TFORM [12]. Since the resulting formulas describing the general second-order stress-energy tensor 
are lengthy and rather complex we shall not display them here. 

The thus obtained approximate stress-energy tensor can be applied in any spacetime provided 
the temporal changes of the geometry are small and A c /L <C 1. The effective action approach that 
we employ in this paper requires the metric to be positively defined. Consequently, the stress- 
energy tensor can be obtained by analytic continuation of its Euclidean counterpart at the final 
stage of calculations. 

The first order (i.e. m~ 2 ) approximation to the renormalized stress-energy tensor of the massive 



scalar, spinor and vector field in the general spacetime has been constructed in Refs 
results generalize the analogous results obtained earlier by Frolov and Zel'nikov 



13 



Jj. This 



161 ] for the 



vacuum type-D metrics as well as the analytic approximation obtained by Anderson, Hiscock and 
Samuel (AHS) for the massive scalar field in a general static and spherically-symmetric geome- 
tries [171 ] . (See also Popov's paper 18].) The AHS approximation is equivalent to the Schwinger 
- DeWitt expansion; to obtain the lowest (i. e. m~ 2 ) terms, one has to use sixth- order WKB 
expansion of the mode functions. 

The range of applicability of such a stress-energy tensor is dictated by the limitations of the 
validity of the renormalized effective action. Numerical calculations reported in Refs. [itI . [l^ ] 
confirm that the Schwinger-DeWitt method provide a good approximation of the renormalized 
stress energy tensor of the massive scalar field with arbitrary curvature coupling as long as the 
mass of the field remains sufficiently large. 
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The stress-energy tensors constructed in Refs. 13|, Il4j. 



physically interesting cases, such as various black holes 
and wormholes 24]. 



171 have been a ppl ied in a number of 



14 



19, 



20 



21 



2211 their interiors 



23] 



In this note we shall calculate the renormalized stress-energy tensor of the massive scalar field 
(in a large mass limit) with an arbitrary curvature coupling in the geometry of the Schwarzschild 
black hole up to m~ 4 terms. We shall also analyze the back reaction problem and briefly study 
the quantum corrected Schwarzschild black hole. Throughout the paper a natural system of units 
is adopted, although in some formulas the constants h, c and G have been, for clarity, restored. 

II. THE STRESS-ENERGY TENSOR 



Now let us return to Eq.([T]) and retain only the first two terms. The approximate stress-energy 
tensor constructed from the coefficients [0,3] and [04] is, therefore, given by 



nab 



1 



1 



5 



(2) 
ab ' 



(5) 



32vr 2 m 2 Jg 5g ab J 32vr 2 m 4 Jg 5g ab 

Since the coefficients [0,3] and [04] are, respectively, the operators of dimension six and eight con- 
structed from the Riemann tensor, its covariant derivatives up to some prescribed order and con- 
tractions, the result of the functional differentiation of the effective action with respect to the 
metric tensor is rather complicated. Moreover, one expects that any attempt to employ the thus 
obtained results for a concrete line element would be, computationally, a real challenge 34|. For 



example, for a general static and spherically-symmetric geometry described by a line element of 
the form 



g tt (r)dt 2 + g rr (r)dr 2 + r 2 {d9 2 + 



sin 



(6) 



the expression describing the next-to-leading term of the stress- energy tensor, when fully expanded, 



consists of 2582 primitive terms for T^ 2t \ 2026 for T^'' and 2066 for T^ ,a . This can be contrasted 
with the number of the primitive terms in the tensor T^ b : 615 for T^ 1 463 for and 634 

Tg^. Fortunately, the final result for a simple metric is, as we shall see, quite simple. 

Making use of the first-order approximation of the stress-energy tensor in the Schwarzschild 
geometry, one easily obtains [if 



(2)r 



,(2)0 



,(1)* 



r (l)r 



327r 2 m 2 r 8 



16 



176 M 

5 r 



V 



M 2 



32ir 2 m 2 r s 



48 M 32 



5 r 



19 626 M 
21 + 315 r 



22 M 1 

45 r + 3 



(7) 



(8) 
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and 



M 2 



32ir 2 m 2 r 8 



22AM 96 \ 734 M 

1 IV H 

5r 5 / 315r 



1 



(9) 



where = £ — 1/6. 

Now, let us consider the second term of the equation ([5]). The second-order calculations are, 
of course, more involved. Fortunately, there are massive simplifications for the Ricci-flat geometry 
and the final result in the Schwarzschild geometry is quite simple. Tedious but routine calculations 
give 



,{2)t 



M 2 



32vr 2 m 4 r 10 



752M 6596M 2S 
144 + j rj 



44 22664M 
Y + 525r 



27166M 2 
525r 2 



(10) 



T (2)r 



M 2 



327r 2 m 4 r 10 



288 1164M 2 1208M\ 12 776M 5506M 2 



7 



7r 2 



+ 



7r / 



r/ + 



75r 



+ 



525r 2 



and 



,(2)0 



M 1 



32vr 2 m 4 r 10 



/ 7760M 2 
V 7r 2 



6084M 1152 \ 1304M 
+ -=- )r/ + 



7r 



25r 



35698M 2 
525r 2 



48 
5 



(11) 



(12) 



The constructed tensor is covariantly conserved, regular and it can easily by checked that at the 



event horizon one has T t 



(2)t 



,(2)r 



Moreover, it should be noted that although the general 



expression describing [a^] involves the terms up to £ 5 , the final result is linear in £. Although, 
generally speaking, there are no limitations placed on the parameter £, two of its values are partic- 
ularly appealing, namely i] = and r] = —1/6 which lead to the conformal and minimal couplings, 
respectively. 

In Figs. 1-3 the run of the (rescaled) components of the stress-energy tensor T b as functions 
of z = r/M for a few exemplary values of the coupling parameter from the range < ^ < 1/6 
is displayed. Although there are strong dependence on £, some general features are common for 

(2)t 

all the curves. Indeed, the T t is negative at the event horizon and remains so for z < 2.1 and 
attains a (positive) maximum. Subsequently it decreases with r approaches a (negative) minimum 

(2)r 

and falls to zero. Inspection of Fig. 2 shows that T r is negative at the event horizon, approaches 
a (positive) maximum and fall to zero as r — > oo. Finally, the run of the angular component (Fig. 
3) is qualitatively similar to that of T t . The behavior of the stress-energy for more exotic values 
of the coupling parameter can easily be inferred form the general formulas (|10til2[) . Specifically, at 
the event horizon one has 

1 



n 



(2)t 



T (2)r 



448007r 2 m 4 (2M) ! 



(12507/ - 53) 



(13) 
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FIG. 1: This graph shows the rcscalcd T^ 1 [A = (8M) 4 7r 2 m 4 ] component of the stress-energy tensor of the 
massive scalar field as a function of z = r/M plotted for a few exemplary values of the coupling parameter 
£. Top to bottom (at the maximum) the curves are plotted for £ = 0.2i (i = 0, 8) and for £ = 1/6. 



and 



T (2)0 _ 

J- a — 



(1500?? - 109) . 



(14) 



268807r 2 m 4 (2M) 8 

Thus far we have carried out our calculations using the Planck units. It is of some interest to 
restore the constants c, G and h in the final expressions describing the renormalized stress-energy 
tensor. Simple manipulations give 

Tjf>" = A {t) x f^ a (z), (15) 

where = G 2 h 3 M 2 /c 5 m 2 r 8 , A( 2 ) = G 2 h 5 M 2 /c 7 m 4 r 10 , and f%\ a {z) are dimensionless functions 
of z = GM/c 2 r. 

Since the Schwinger-DeWitt approximation is local and the geometry at the event horizon is 
regular, one expects that the stress-energy tensor is also regular there. On the other hand, the 
stress-energy tensor is regular in the physical sense if it is regular in a coordinate system which is 
well behaved as r — > r + . For example, the components of the stress-energy tensor in a freely 
falling frame, denoted here as T( )( ), 7"(o)(i) an d ^(i)(i) are 

7 2 (t;-t/) 



(0)(0) 



■(1)(1) 



/ 

f 



- T r 

r i 



+ T r 



(16) 
(17) 
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FIG. 2: This graph shows the rescaled Tr [A = (8M) 4 7r 2 m 4 ] component of the stress-energy tensor of the 
massive scalar field as a function of z = r/M plotted for a few exemplary values of the coupling parameter 
£. Top to bottom (at the maximum) the curves are plotted for £ = 0.2i (i = 0, 8) and for £ = 1/6. 



(0)(1) 



7Vt^7(t;-t/) 



(18) 



where 7 is the energy per unit mass along the geodesic and f(r) = —gtt( r )- Inspection of Eqs. (1161 
TBj) shows that if all components of T% and (TJT — T*)/ / are finite on the horizon the stress-energy 
tensor in a freely falling frame is finite as well. 

Now, simple calculations show that the difference between radial and time components of the 
stress-energy factors 



where 



T Ut_ T p r= { 1 _^\ F (i) ij 



(19) 



M 1 



7r 2 m 2 r 8 \ 10 



7 



13 

336 



(20) 



F (2) (r) 



M 2 



7r 2 m 4 r 10 



81 
14 



485 M 



7 1021 M 
20 + 1050 T 



(21) 



28 r 

i = 1,2, and, consequently, both tensors are regular in a physical sense. Moreover, using our 
general formula describing the stress- energy tensor it can be shown that it remains so in any static 
and spherically-symmetric spacetime. 
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FIG. 3: This graph shows the rescaled T g [A = (4M) 8 7r 2 m 4 ] component of the stress-energy tensor of the 
massive scalar field as a function of z = r/M plotted for a few exemplary values of the coupling parameter 
£. Top to bottom (at the maximum) the curves are plotted for £ = 0.2z (i = 0, 8) and for £ = 1/6. 



Having constructed the next-to-leading term of the renormalized stress-energy tensor which 
depends on a general metric one can analyze the back reaction of the quantized field upon the 
black hole geometry. It should be emphasized once more that accepting the approximation ([SJ we 
ignore particle creation which is a nonlocal effect. To simplify our discussion we shall assume that 
the cosmological constant and the renormalized coupling parameters a and f3 in the quadratic part 
of the total action 



where C^cd is the Weyl tensor, identically vanish. The semiclassical Einstein field equations have, 
therefore, a standard form 



where (T a b[g]) = 0(h) is the renormalized stress-energy tensor. 

Since the total stress-energy tensor depends functionally on a wide class of metrics, one can, in 
principle, construct the self-consistent solution of the system (|23p . It should be noted however, that 
since the general stress-energy tensor, T a , is constructed from [04] it contains the terms up to eight 



III. THE BACK REACTION PROBLEM 




(22) 



G b a [g]=8n{T ab [g]), 



(23) 



9 



derivatives of g a b, and, consequently, there is a real danger that the semiclassical equations may lead 
to physically unacceptable solutions [25]. Moreover, the tensor (T%) is extremely complicated and 
it is natural that one is forced to refer to some approximations. Here we shall treat the right hand 
side of the semiclassical Einstein field equations as perturbation. Restricting to the perturbative 
solutions of the effective theory may be, therefore, the only one way to obtain the (approximate) 
physical solutions. 



For the quantized massless fields in the Schwarzschild geometry the back reaction program 



been initiated by York 



221, 



28 



has 



29 



30, 



261 ] . Subsequently, it has been applied in numerous papers 
31], where various aspects of the back reaction of the massive fields upon the black hole geometry 
has been studied using the first order approximation to the stress-energy tensor. 

In the semiclassical approach we ignore the effects caused by the quantized gravitational field 
simply because they are not known. On general grounds, however, it is expected that the pertur- 
bations of the classical metrics caused by gravitons should be of the same order as from the other 
quantized fields. One can justify neglecting of the graviton contribution to the total stress-energy 
tensor taking a large number of various fields in the calculations. It can be argued that the contri- 
bution of the gravitons would be small as compared with the total effect caused by other physical 
fields. 

n 

Now, let us introduce the dimensionless parameter e [32] and make the substitution (T a b[g]) — > 
e(T a b[g]). Expanding the metric tensor as 

9a b = gf b ] +eg^+0(e% (24) 



inserting it into the semiclassical equations (|23p and collecting the terms with the like powers of 
the auxiliary parameter, one obtains 

G b a [g^} = (25) 

and 

G^}=8,(tWV^+tP%^]), (26) 



i.e., the modifications of the geometry caused by stress-energy tensor calculated in the corrected 
black hole spacetime, T^ b [g^], are ignored as these additional terms would be 0(h 2 ). In other 
words we are looking for 0(H) corrections to the classical solution. 

From Eqs. ([Ti ll 2 j) one sees that the solution of the back reaction problem reduces to elementary 
quadratures. First, let us consider the issue of the integration constants which appear in solutions 
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of the differential equations. The zeroth-order equations will yield two integration constants, say, 
ci and C2, which can be set to —1 and M, respectively. Here M is a "bare" mass of the black hole 
and ci can be determined from the condition g$ 9rr = — 1- On the other hand, the integration 
constant C\ appearing in the component g$ of the metric tensor can be absorbed in a process of 
the finite renormalization of mass. Indeed, it can be demonstrated that with the substitution 

M = M - ^ed (27) 

the radial component of the metric tensor can be written as 

97r l (r) = 1 - m + *>«(r) + 7f )(r) + 0(e 2 ), (28) 

where \r) and Vr 2 \r) are given respectively by 

, {22M \ 313M 19 , , 



and 



2100r 35 9450r 2 V 7 63r2 2r / 

Similarly, for g a to 0{h) one has 



, 2833A4 11 135837W 2 /36 16497W 2 47A4 \ , , 

^ 2) (r) = -^Tr^r + « + " hr + -55=5 sr- W ( 30 ) 



<teM = -^(r) exp (2e^(r)) , (31) 



where 



^ ~ vrm 2 r 6 V 15 77 ~ 504/ + vrm 4 r 8 



2042A4 7 /4857W 81 
4725r ~ 40 ~ V 63r 28 



+ C 2 . (32) 



The integration constant C 2 can be fixed by demanding that the time component of the metric 
tensor approaches its Minkowskian value as r — > 00, that is equivalent to normalizing the time 
coordinate at infinity. 

As before, it is of some interest to restore the physical constants. Putting M = GAi/c 2 , Ipi = 
(h/Gc 3 ) 1 / 2 and A c = h/mc in Eq. (|28p the radial component of the metric tensor can schematically 
be written as 

, M 2 X 2 l 2 , x M 2 \tl 2 , , s 
g- 1 = l-2z + e -±PLWi{z;r))+e ^-^ W 2 {z;r,), (33) 

where W\ and W 2 are simple polynomials depending parametrically on 77 and their exact form can 
easily be inferred from Eq. (j33"j) . A similar expression can be constructed for g u , and the result 
can be schematically written in the form 

j,6 



g tt = -l + 2z + e - \ c Pl V 1 (z; V ) + e ^ V 2 (z; V ), (34) 
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where V{ comprise another pair of simple polynomials. 

The location of the event horizon of the quantum-corrected Schwarzschild black hole is de- 
termined by the equation g tt (r + ) = 0. Putting r + = r+ + er^ 1 one concludes that r + is given 
by 

r + = 2M H ^i- 77 H c 77 . (35) 

+ 480vrM 3 V 504 J 2016vrM 5 V1200 'J v ' 

The Euclidean version of the line element ([6]) obtained with the aid of the Wick rotation has no 

conical singularity provided the "time" coordinate is periodic with a period (5 given by 

= ftm + 4vr {g tt grr) 1/2 (Jj~9tt) ■ (36) 
For the quantum-corrected metric (|28H32|) . the period is 



^ 8irM 302407W 3 m 2 1512007W 5 m 4 (3?) 



to the first order in e. 

The surface gravity, k, which is proportional to the temperature of the black hole can be 
calculated (for the Lorentzian metric) from a simple relation 

£ |r=r_)_ 

where k a is a timelike Killing vector. The Hawking temperature, Th = k/2it, is, therefore, given 
by 



^(-9ug rr ) |— | )r=r+ - - 



1 ( l 1 1 



8ttM TT 2 m 2 (AM) 5 \1980 9450.M 2 m 2 
It should be noted that the temperature, Tjj, when expressed in terms of the total mass of the 
system as seen by a distant observer, is independent of the coupling constant. 

On the other hand, one can express the results (|28ti32|) in terms of the horizon defined mass, 
Mjj = r + /2 which, of course, differs from the total mass of the system as seen by a distant observer. 
It can be achieved, for example, by inverting Eq. (|35p and the elementary manipulations give 



M = M H ( v - —) £ —r-^ ( — - ryV (40) 

9607rm 2 M| V 504 / 40327rm 4 M| V 1200 'J v ' 

Consequently, one can systematically substitute Eq. (f4T)|) in Eqs. (|28M32|) . expand and finally linearize 
the thus obtained results. 
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Equally well, one can start with a slightly different representation of the line element putting 

9 „ W = -exp(2^))(l-^V = (41) 
\ r J r 

expanding the functions M(r) and ip(r) into the power series 

k 

M(r) = ^M i {r)e i (42) 

k 

^{r)=Y J ^ i (43) 

i=l 

and retaining only the linear terms. Now, accepting the boundary condition ^i(oo) = 0, repeating 
the calculations for the line element (|41|) . with the integration constants appearing in the solution 
for M(r) determined either from Mq(oo) = M and Mi(oo) = or from Mo(r + ) = r + /2 and 
Mi(r+) = 0, one can easily reconstruct all our previous results. 



IV. FINAL REMARKS 



In this paper we report our calculations of the next-to-leading term of the renormalized stress- 
energy tensor of the quantized massive scalar field in a large mass limit. To achieve this, we have 
calculated the effective action constructed form the (integrated) coincidence limit of the coefficients 
as(x,x') and (24(2, x'), and, subsequently, we have calculated the approximate stress-energy tensor 
by functional differentiation of the thus obtained action with respect to the metric tensor. The 
obtained stress-energy tensor can be employed in any spacetime provided the condition Xc/L <C 1 
holds. The general formulas describing the stress-energy tensor are extremely complex, but, when 
applied to the Schwarzschild geometry, they yield remarkably simple result, which is the main 
result of this paper. The general stress-energy tensor has been used in the analysis of the back 
reaction of the quantized field upon the geometry of the Schwarzschild black hole. 

Finally, we indicate a few possible directions of investigations. First, it would be interesting 
to examine the vacuum polarization effects in more complex backgrounds, as for example, the 
spacetime of the electrically charged black holes with or without the cosmological constant. Further, 
the numeric approach to the back reaction would certainly strengthen our understanding of the 
problem. This group of problems is actively investigated and the results will be published elsewhere. 
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